Bell's theorem for general N-qubit states 
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We derive a single general Bell inequality which is a necessary and sufficient condition for the 
correlation function for N particles to be describable in a local and realistic picture, for the case in 
which measurements on each particle can be chosen between two arbitrary dichotomic observables. 
We also derive a necessary and sufficient condition for an arbitrary N-qubit mixed state to violate this 
inequality. This condition is a generalization and reformulation of the Horodeccy family condition 
for two qubits. 
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Local realism imposes constraints on statistical corre- 
lations of measurements on multiparticle systems. They 
are in the form of Bell-type inequalities ||, ^, ^, ||, 
^, I?! ||. In a realistic theory the measurement results 
are determined by "hidden" properties the particles carry 
prior to and independent of observation. In a local real- 
istic theory the results obtained at one location are in- 
dependent of any measurements, or actions, performed 
at space-like separation. Quantum mechanics predicts 
violation of these constraints. This is known as Bell's 
theorem [|] . 

However the problems a) what are the most general 
constraints on correlations imposed by local realism, and 
b) which quantum states violate these constraints, are 
still open. The latter has been solved in general only in 
the case of two particles in pure states 10 1 and for 
two-qubit mixed states Only recently bounds for 

local realistic description of a higher-dimensional system 
have been found in some simple cases [|l2|, |l3|, 

Here the answer to the two long-standing questions 
(a) and (b) is presented for the case of a standard Bell 
type experiment on N qubits. By a standard Bell exper- 
iment we mean, one in which each local observer is given 
a choice between two dichotomic observables. We first 
derive a single general Bell inequality that summarizes 
all possible local realistic constraints on the correlation 
functions for a N-particle system. From this inequality 
one obtains as corollaries the Clauser-Horne-Shimony- 
Holt (CHSH) inequality Q for two-particle systems 
and the Mermin-Ardehali-Belinskii-Klyshko (MABK) in- 
equalities for N particles [|[ ||, |^. We show that the 
correlation functions in a standard Bell experiment can 
be described by a local realistic model if and only if the 
general Bell inequality is satisfied. Therefore the general 
Bell inequality is a sufficient and necessary condition for 
correlation functions, in such an experiment, to be de- 
scribable within a local realistic model. We also find a 
necessary and sufficient condition for correlation func- 
tions for N qubits in an arbitrary (mixed) quantum state 
to violate the general Bell inequality in direct measure- 
ments. This condition is generalization and reformula- 
tion of the one given by the Horodeccy family llll for 



two qubits. 

These results are not only of importance from the fun- 
damental point of view, but also as a research towards 
identifying ultimate resources for quantum information 
processing. Recently it was shown that there is a 
direct link between the security of the quantum commu- 
nication protocols, and the violation of Bell inequalities. 

We shall now derive the general Bell inequality. Con- 
sider N observers and allow each of them to choose be- 
tween two dichotomic observables, determined by some 
local parameters denoted here ni and n2- We choose 
such a notation of the parameters for brevity; of course 
each observer can choose independently two arbitrary di- 
rections. The assumption of local realism implies ex- 
istence of two numbers Aj{ni) and Aj(n2) each taking 
values -f-l or -1, which describe the predetermined result 
of a measurement by the j-ih. observer of the observable 
defined by fti and n2, respectively (we do not discuss 
stochastic hidden variable models, as they always can be 
constructed from underlying deterministic ones). In a 
specific run of the experiment the correlations between 
all N observations can be represented by the product 
W^^i Ajiji};.), with kj = 1,2. The correlation function, 
in the case of a local realistic theory, is then the average 
over many runs of the experiment 
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The following algebraic identity holds for the predeter- 
mined results: 
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J2 Sisi, ...,SN)Y[[Aj{ni) + sjAj{ri2)] =±2 

si,..,sjv = ±l j=l 

(2) 

where S{si, sn) stands for an arbitrary function of 
the summation indices si, sn G {—1, 1}, such that its 
values are only ±1, i.e. S{si, sn) = ±1. To prove this 
identity, note, that, since Aj{n) = ±1, for each observer 
j one has either \Aj{ni) + Aj{n2)\ = and \Aj{ni) — 
Aj{n2)\ = 2, or the other way around. Therefore, for all 
sign sequences of si,...,sn the product njLi[^i("-i) + 



2 



SjAj{n2)] vanishes except for just one sign sequence, for 
which it is ±2^. If one adds up all such 2^ products, 
with an arbitrary sign in front of each of them, the sum is 
always equal to the value of the only non- vanishing term, 
i.e., it is ±2^. 

After averaging the expression (||) over the ensemble of 
the runs of the experiment (compare Eq. (|l|)) one obtains 
the following set of Bell inequalities 

I J2 S{s,,...,Sn)J2 S^^-^.4"~^i?(fcl,...,fcAr)| <2^. 

= = « fci,...,fcjv 

= -1,1 —1,2 

(3) 

Since there are 2^ different functions S{si, sn), 
the inequalities (^) represent a set of 2^ Bell inequal- 
ities for the correlation functions. Many of these in- 
equalities are trivial. E.g., when the choice for the func- 
tion is S'(si, sjv) = 1 for all arguments, we get the 
condition £^(1,1, ...1) < 1. Specific other choices give 
non-trivial inequalities. For example, for S{si, sn) = 
V^cos{-j + {si + ... + sn-N)j) one recovers the MABK 
inequalities, in the form derived by Belinskii and Klyshko 
Specifically, for N — 2, the well known CHSH in- 
equality follows. For = 3, one obtains the inequal- 
ity 

\E{1, 2, 2) -I- E{2, 1, 2) + E{2, 2, 1) - E{1, 1, 1)| < 2. (4) 

Inequalities, like the one above, with the minus sign at 
a different location, and/or measurements 1 and 2 per- 
muted, form together an equivalence class. 

The full set of all 2^ inequalities is equivalent to 
the single general Bell inequality ^ |l8| 

J2 I E Eih,...,kN)\<2''. (5) 

= -1,1 =1,2 

The equivalence of (||) and (j^) is evident, once one recalls 
that, for real numbers one has |a-|-6| < c and |a — 6| < c if 
and only if |a| -I- 16| < c, and writes down a generalization 
of this property to sequences of an arbitrary length. 

Thus far we have shown that when a local realis- 
tic model exists, the general Bell inequality (|^) follows. 
The converse is also true: whenever inequality (^) holds 
one can construct a local realistic model for the correla- 
tion function, in the case of a standard Bell experiment. 
This establishes the general Bell inequality (||) presented 
above as a necessary and sufficient condition for local 
realistic description of A^ particle correlation functions 
in standard Bell-type experiments. This is why one can 
claim that the set of Bell inequalities (||) is complete. 

The proof of the sufficiency of condition (^) will be 
done in a constructive way. A local realistic theory must 
ascribe certain probabilities to every possible set of pre- 
determined local results. Just like if the local measuring 
stations were receiving instructions, what should be the 
measurement results for (here) two possible settings of 
the local apparatus. 



One can ascribe to the set of predetermined local re- 
sults, which satisfy the following conditions Aj{ni) = 
SjAj{n2), the hidden probability 

p{si,..,SN) = ^\Y.'l""'-'^N"''E{h,...,kN)\, (6) 
fci,..,fcjv 

and one can demand that the product Ilj^^i ^ji^i) 
the same sign as that of the expression inside of the mod- 
ulus defining the p{si, sn). In this way every definite 
set of local realistic values is ascribed a unique global hid- 
den probability. However, if the inequality (^) is not sat- 
urated the probabilities add up to less than 1. In such a 
case, the "missing" probability is ascribed to an arbitrary 
model of local realistic noise (e.g., for which all possible 
products of local results enter with equal weights). The 
overall contribution of such a noise term to the correla- 
tion function is nil. In this way we obtain a local realistic 
model of a certain correlation function. 

However, one should check whether this construction 
indeed produces the model for the correlation function 
for the set of settings that enter inequality (||), that is 
for E{ki, kN)- For simplicity take N ^ 2. One can 
build a "vector" {E{1,1), E{1,2), E{2,1), E{2,2)) out of 
the values of the correlation function. The expansion co- 
efficients of this "vector" in terms of the four orthogonal 
basis vectors (1, si, S2, S1S2) (recall, that si, S2 € {—1,1}) 
are equal to the expressions within the moduli entering 
inequality (||). By the construction shown above the local 
realistic model for A^ = 2 gives the following "vector" 

(i?Lfl(l, 1), i?Lfl(l, 2), i?Lfl(2, 1),Elr{2, 2)) (7) 

= i E [ E sJ^-^4^-^i?(fcl,fe)](l,Sl,S2,5lS2). 
Sl,S2 ki,k2 

Since the vector built out of the correlation function val- 
ues and its local realistic counterpart have the same ex- 
pansion coefficients in the basis, they are equal. Thus, 
the sufficiency of (|5]) as a condition for the existence of a 
local realistic model is proven. The generalization to an 
arbitrary A^ is obvious. 

Quantum mechanical predictions can violate the in- 
equality (H). Simply, if a MABK inequality is vio- 
lated, then the general inequality, which also includes the 
MABK inequalities, is violated too. However, the con- 
verse statement is not always true. The new inequality is 
more restrictive. In the problem of identifying quantum 
states of highly nonclassical traits, it is important to find 
the class of quantum states, which are not describable by 
local realistic models. We will now derive the necessary 
and sufficient condition for an arbitrary (pure or mixed) 
quantum state to violate the general Bell inequality (^. 

An arbitrary mixed state of A^ qubits can be written 
down as 

1 ^ 

E T^i--N '^li®-®'^^^, (8) 

X\,...,XJS[—^ 
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where Ug is the identity operator in the Hilbert space of 
qubit j , and a^^ are the Pauh operators for three orthog- 
onal directions Xj = 1,2,3. The set of real coefficients 
Txi...xNJ with Xj — 1,2,3 forms the so-called correlation 
tensor T. The correlation tensor fully defines the N-qubit 
correlation function: 

EQM{ki,...,kN) = Tr[p {fik^ ■ a (g) ...(g) fikj^ ■ a)] (9) 

3 

(10) 

a; 1 , . . . , a: 71 — 1 

where {nkj)xj are the three Cartesian components of the 
vector ftk^ . For convenience we shall write down the last 
expression (^0|) in a more compact way as (T, n^.^ ... (g) 
fikfi), where (..., ...) denotes the scalar product in R^^ . 

We now insert the quantum correlation function 
EqAiiki, fcAr) into the Bell inequality (^, and obtain 

2 2 

(11) 

This inequality can be simplified. For each observer there 
always exist two mutually orthogonal unit vectors d\ and 
and the angle Uj such that X]fcj=i ^fcj ~ '^^i cos(aj + 
^) and X]fcj=i(~l)''^^fc3 — 2a^cos(Q!j + tt). Using the 
notation c^. —cos{aj +Xj^), one can write the inequality 
(P) as 



J2 |4,...cf„(f,ai^^...^af„)|<l. (12) 
One can transform this inequality into 

2 



(13) 



Xi,...,Xn — ^ 



where Tx-^...xn is now a component of the tensor T in a 
new set of local coordinate systems, which among their 
basis vectors have a{ and . The two vectors serve as the 
unit vectors which define, say, the local directions x and 
y. The values of cl_ enter (|l^) directly, not as moduli, 
because without this constraint the maximal value of the 
left hand side does not change. 

We conclude from the above reasoning that the nec- 
essary and sufficient condition for an arbitrary N-qubit 
state to satisfy the general Bell inequality ^ can be put 
in the following way. The correlations between the mea- 
surements on N qubits satisfy inequality if and only 
if in any set of local coordinate system of N observers, 
and for any set of unit vectors c? = (c{,c^) one has 



mod 



E 

X\,...,X j^ — l 



IT 

'-xx---'^xn 



< 1. 



(14) 



Let us give a geometric interpretation of ([T^). Sup- 
pose one replaces the components of the correlation ten- 
sor Txi...xN by their moduli \Txi...xn\^ and builds of such 
moduli a new tensor T"^°'^. Suppose moreover one trans- 
forms this modified tensor into a new set of local coordi- 
nate systems, each of which is obtained from the old one 
by a rotation within the plane spanned by axes 1 and 
2 of the initial coordinates. If this new tensor satisfies 
constraint ([l^) for an arbitrary choice of the initial set 
of local coordinate systems, then, and only then, a local 
realistic description of correlation function is possible, in 
the case of any standard Bell experiment. 

In other words, T™""^^^ is a component of T"''°'^ along 
directions defined by the unit vectors cP , j = 1, N . If 
the condition (^4|) holds, then the transformed compo- 
nents T^°'^^^ do not have values larger than 1. Only 
then, they can describe products of local results, which 
are only of the values ±1, like is for any correlation ten- 
sor. One therefore can express the condition ([ij) as fol- 
lows: within local realistic description T™°'' is also a pos- 
sible correlation tensor. This bears a similarity with the 
Peres [T9[ | necessary condition for separability (a partially 
transposed density matrix is a possible density matrix). 
Note that (|lj) could also be put in yet another way: ar- 
bitrary changes of the signs of some of the coordinates of 
T still leave it as a possible correlation tensor. 

By applying the Cauchy inequality to the middle term 
of expression ( P ) one obtains directly the following use- 
ful and simple sufficient condition for local realistic de- 
scription of the correlation functions for N qubits. If in 
any set of local coordinate systems of N observers 



E 



< 1 



(15) 



then the correlations between the measurements on N 
qubits satisfy the general inequality (H). 

By performing rotations in the planes defined by direc- 
tions 1 and 2 of each of the N observers one can vary the 
values of the elements of the correlation tensor, but these 
variations do not change the left-hand side of inequality 
(|l|). In this way, one can find local coordinate systems 
for which some of the correlation tensor elements vanish. 
Thus the criterion ( |l5| ) can involve a smaller number of 
them (compare the three qubit case in Re f M ). 

There are special situations for which (|l5) turns out 
to be both the necessary and sufficient condition for cor- 
relation functions to satisfy the general Bell inequality 
(h). Formally this arises whenever the two "vectors" 



{cl 



.c^) and (|Tn...i|,...,|r22...2|) in (|T|) are 



parallel. Only then, since the first vector has a unit norm, 
the expression on the left hand side of (|l^) reaches the 
one on the left hand side of (^5|) and thus the conditions 
( P ) and (|l^) are equivalent ones. 

Let us consider two examples of application of our 
results. We first study an arbitrary two-qubit state 
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to recover the Horodeccy condition In this 

case, since the two "vectors" (|Tii|, jTial, |T2i|, IT22I) and 



{c\c\,c\c2,c\c\,c\c\) in (|T^) can be made parallel by a 
suitable choice of free parameters, (|l^) is the necessary 
and sufficient condition for the violation of local realism. 
Two of the parameters come from the arbitrariness in se- 
lecting the two local coordinate systems (i.e. they come 
from arbitrary transformation of the correlation tensor 
by rotations within 1-2 planes of each of the two local 
coordinate systems). Two more parameters are the two 
angles aj which define the second "vector". Therefore 
the condition reads Tf^ + T^^ + T^i + T^2 < 1- In ad- 
dition one can always find local coordinate systems such 
that T12 = T21 = and our condition transforms into 
+ < 1 ■ This is equivalent to the Horodeccy con- 
dition , as T^^ and for the diagonalized correlation 
tensor are equal to two eigenvalues of the matrix T"^T, 
where is the transposed T. 

As another example, we consider the Werner states. 
Such states have the form pw = V\4'ghz){''Pghz\ + {^ — 
V)pno^se, where \iPghz) = ^{\Q)i-\Q)n + \1)i-\1)n) is 
the maximally entangled (GHZ) state [|| and pnoise ~ 
1/2^ is the completely mixed state. Here the weight V 
of the GHZ-state can operationally be interpreted as the 
interferometric contrast observed in a multi-particle cor- 
relation experiment. The nonvanishing components of 
the correlation tensor in the xy planes for the Werner 
state are: the components which contain an even num- 
ber of y's and Txx...x- There are alltogether 2^""'^ of 
them. Their values are either -\-V or ~V . Since again 
the two "vectors" in (O) can be made parallel, (15) is 



the necessary and sufficient condition for the violation 
of local realism. Indeed, if one rotates all but one local 
coordinate systems by 45°, then all 2^ components of 
the "vector" (|rii...i|, IT22...2I) become equal to V/V2. 
Furthermore, if one chooses all aj equal to — 7r/4, the 
unit "vector" {c\...Ci , ...^c\...C2) has all its components 
equal to 1/V2^. Therefore, the two "vectors" are paral- 
lel. Thus, using criterion (|5|) we conclude that the corre- 

latirjn fiinrtinns for ihp Wprnpr utatp dpfinitnly rannnt bp 



ment choices for each observer than two. One can expect 
in such cases even stronger restrictions for the local real- 
istic description 
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dpsrlrihprl by Inral rpalism if and only if V > 1 /^/9^~l 



More applications of the formalism, leading to some 
unexpected results, are given in [pTf . There a family of 
pure entangled states is found, which do not violate any 
Bell inequality for correlation functions, for the standard 
Bell experiment. 

It will be interesting to see generalizations of the cri- 
teria for violation of local realism to the cases of higher- 
dimensional systems than qubits and to more measure- 
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